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摘要 : 消除 综合 孔径 望远镜 子 孔径 之 间 的 相对 光 程 差 是 实现 高 分 辩 率 干涉 成 像 的 前 提 条 
件 ， 条 纹 检 测 法 是 一 种 检测 相对 光 程 差 的 有 效 办 法 。 因 子 孔 径 的 空间 位 置 排 布 使 干涉 条 纹 具 
有 一 定 的 方向 性 ， 若 不 能 精确 获知 干涉 条 纹 的 角度 ， 则 无 法 沿 条 纹 的 法 线 方向 行 采样 ， 进 而 
无 法 根据 对 比 度 变 化 曲线 的 最 大 值 获得 子 孔径 之 间 的 最 小 光 程 差 位 置 。 提 出 了 一 种 基于 傅 里 
叶 变 换 和 图 像 二 值 化 阅 值 遍历 的 获得 干涉 条 纹 角 度 的 方法 ， 首 先 介绍 了 算法 的 基本 原理 ， 其 
次 通过 对 条 级 角 度 为 43 "的 仿真 数据 进行 算法 验证 得 到 的 角度 为 43.007 8"， 与 理论 值 的 误 
差 为 0.018% ， 证 实 了 方法 的 可 行 性 。 最 后 对 比 了 未 旋转 相机 和 旋转 相机 两 种 情况 下 的 条 纹 
对 比 度 变 化 曲线 ， 可 知 通 过 旋转 相机 使 条 纹 转 至 相机 靶 面 纵 轴 方 向 再 进行 采样 的 办 法 ， 更 有 
利于 精确 得 到 相对 光 程 差 的 最 小 位 置 。 
关键 词 : 傅 里 叶 变 换 ; 阅 值 ; 干涉 条 纹 ; 条 纹 角 度 ; 光 程 差 
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对 空间 目标 高 分 状 率 成 像 一 直 以 来 都 是 天 文 工 作者 不 断 努力 的 目标 ， 在 波长 一 定 的 条 件 下 ， 望 远 
镜 的 分 辨 率 通常 受 限于 系统 的 孔径 大 小 ， 孔 径 越 大 分 辨 率 越 高 ”， 但 材料 、 工 艺 及 成 本 的 因素 限制 
了 系统 孔径 的 增 大 “”” 。 综 合 孔径 成 像 技 术 可 以 通过 容易 加 工 的 小 孔径 合成 大 孔径 系统 ， 实 现 高 分 辨 
率 成 像 ， 但 前 提 条 件 是 消除 子 孔径 间 的 相对 光 程 差 。 条 纹 检 测 法 是 检测 综合 孔径 望远镜 子 孔 径 之 间 相 
对 光 程 差 的 一 种 有 效 方法 “。 通 常 在 利用 相 移 法 得 到 综合 孔径 两 孔径 之 间 的 干涉 条 纹 后 ， 需 要 对 
连续 采集 的 多 组 条 纹 数 据 行 采样 ， 从 而 得 到 最 小 光 程 差 的 位 置 。 但 因为 子 孔径 的 空间 排 布 并 不 是 上 下 
或 左右 并 列 排 布 ， 因 而 得 到 的 条 纹 数据 有 一 定 角度 。 大 得 到 条 纹 的 角度 ， 则 可 以 沿 干 涉 条 纹 的 法 线 方 
向 采样 ， 获 得 对 比 度 的 变化 曲线 ， 从 而 根据 对 比 度 的 峰值 位 置 找到 子 孔径 之 间 的 最 小 相对 光 程 差 位 
置 ， 最 终 实现 综合 孔径 望远镜 的 高 分 辨 率 成 像 。 


1 DEWR 


将 独立 三 孔径 干涉 条 纹 作为 仿真 对 象 ，3 个 子 孔径 的 空间 位 置 呈 等 边 三 角形 排列 ， 如 图 (a), 
理论 上 ， 透 射 到 光束 合成 镜 和 人 瞳 面 上 的 光束 相对 于 透射 到 整个 系统 入 瞳 面 上 的 光束 是 等 比例 复制 的 ， 
但 实际 上 由 于 剪 切 误 差 和 放大 率 误 差 的 原因 ， 两 者 之 间 存 在 瞳 映 射 误 差 ， 所 以 映射 到 合成 镜 和 人 腾 面 上 
的 3 个 光斑 不 是 等 边 三 角形 分 布 ， 其 相对 位 置 如 图 1(b) 。 

由 图 1(b) 可 以 得 到 基线 B1, 与 水 平方 向 的 夹 角 近似 为 -47 (这 里 的 负 号 表示 与 水 平 x 轴 负 方向 的 夹 
fi; 大 为 正 ， 则 是 与 正方 向 的 夹 角 )， 由 于 被 基线 切割 ， 条 纹 与 水 平方 向 的 夹 角 约 为 43”。 当 单口 径 望 
远 镜 对 无 穷 远 处 的 点 目标 成 像 时 ， 在 系统 的 焦 平 面 上 接收 到 的 是 点 目标 的 夫 环 禾 费 衍射 像 ， 因 此 子 孔 径 
1 和 子 孔 径 3 对 无 穷 远 处 的 点 目标 成 像 时 ， 在 系统 的 焦 平 面 上 接收 到 的 衍射 光 强 可 以 分 别 表示 为 
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双 和 孔径 干涉 实质 上 是 在 单 孔径 衍射 的 基础 上 铸 加 了 一 个 与 基线 有 关 的 余弦 调制 项 ， 由 此 可 以 得 到 
双 孔 径 干 涉 条 纹 的 计算 公式 : 


N+AA/2 
2 
I, = | Vly + + 21, x x cos [x Cp x Ba +8) 114A, (3) 
A0-AA/2 


FP, J, AB ZB SER eR x 为 焦 面 位 置 坐 标 ; A, 为 光源 中 心 波长 ，AA 为 光源 带宽 ; f 为 合 
RAEE; 8 为 泊 松 误差 。 根 据 图 1(b) 可 以 计算 3 条 基线 长 度 的 平均 值 有 约 为 16mm， 具 体 仿真 参数 
如 下 : A6=600 nm，AA=100 nm, f=1260mm， 像 元 尺寸 Pixel=7.4 hm， 光斑 大 小 D=9mm，6=0， 仿 
真 条 纹 如 图 2。 


TA A: (0.00, 10.00) 
B: (-5.00, -5.00) 
C: (11.00, -2.00) 
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图 1 (a) 子 孔径 空间 排 布 图 ; (b) 合成 镜 入 瞳 面 光束 相对 位 置 图 
Fig.1 (a) Layout of sub-apertures; (b) Beam position at the entrance pupil of lens 


2 算法 原理 及 仿真 条 纹 计算 


2.1 Hough 变换 直线 提取 法 计算 条 纹 角度 

Hough 变换 直线 提取 是 一 种 比较 经 典 的 算法 | 采用 证 据 收 集 的 方式 ， 遍 历 一 幅 图 像 上 的 所 有 直 
线 的 位 置 ， 哪 条 直线 上 的 特征 点 最 多 ， 哪 条 直线 就 是 期 望 得 到 的 直线 。 利 用 此 算法 思想 计算 条 纹 角 
度 的 主要 步骤 是 : 先 将 条 纹 的 灰 度 图 像 进行 边缘 提取 ， 然 后 检测 边缘 提取 图 像 的 直线 段 ， 根 据 提 取 
的 直线 得 到 条 纹 的 斜率， 从 而 得 到 条 纹 的 角度 。 对 图 2 的 仿真 条 纹 图 进行 Hough 变换 直线 提取 的 结 
如 图 3。 

由 于 条 纹 的 边缘 并 不 完全 是 直线 ， 对 于 比较 细 长 的 条 纹 ， 此 方法 尚 可 ， 寿 条 纹 比 较 宽 、 比 较 短 ， 
此 方法 的 计算 结果 有 很 大 偏差 。 此 外 ， 算 法 在 边缘 检测 和 直线 提取 时 都 需要 用 到 浆 值 ， 浆 值 的 选择 对 
计算 结果 有 很 大 影响 。 传 里 时 变 换 计 算 条 纹 角度 的 方法 不 依赖 于 条 纹 的 长 度 及 宽度 ， 同 时 只 需 在 二 值 
化 时 用 到 阔 值 ， 大 大 简化 了 算法 对 国 值 的 依赖 程度 。 
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T) 边缘 检测 sue 


图 2 仿真 条 级 图 图 3 对 仿真 条 纹 图 的 堆 夫 变换 直线 提取 


Fig.2 Simulation fringe Fig.3 Hough transform line extraction from simulation fringe 


2.2 傅 里 叶 变 换 计算 条 纹 角 度 

对 于 一 个 光学 成 像 系 统 ， 通 常用 系统 的 点 扩散 函 
数 (Point Spread Function, PSF) 描述 点 目标 在 空间 域 
的 成 像 性 能 ; 在 频率 域 则 用 系统 的 光学 传递 函数 
(Optical Transfer Function, OTF) 描述 中 。 健 里 叶 变 
换 汪 是 将 图 像 从 空间 域 转换 到 频率 域 ， 建 立 点 扩散 
函数 与 光学 传递 函数 之 间 的 联系 。 光 学 传递 函数 的 模 
PRA VA till (& 6 PK BL ( Modulation Transfer Function, 


RMA PPRLES 


MTF) ， 整 个 系统 的 调制 传递 函数 由 子 调制 传递 函数 检测 连通 区 域 ， 并 求 每 个 连通 区 域 质 心 
组 成 ， 子 调制 传递 函数 在 频率 域 分 布 的 角度 方向 由 子 
孔径 的 排 布 决定 ， 这 意味 着 如 果 知 道 调 制 传递 函数 在 计算 质心 连 线 的 斜率 及 角度 


频率 域 的 位 置 分 布 ， 就 能 得 到 子 孔径 的 空间 位 置 排 
A?) ， 从 而 得 到 基线 方向 以 及 条 纹 的 角度 。 基 于 此 
理论 ， 本 文采 用 基于 传 里 叶 变 换 求 条 纹 角 度 的 具体 算 
法 见 图 4。 

把 条 纹 图 像 的 灰 度 数据 通过 傅 里 叶 正 变换 转换 成 
图 像 的 频率 分 布 ， 如 图 5(a) 。 根 据 传 里 叶 变换 的 性 
质 ， 原 来 集中 在 图 片 中 心 的 能 量 通过 傅 里 叶 变换 后 ， 
能 量 将 分 布 于 图 片 的 4 个 角 ， 所 以 需要 将 分 布 于 4 个 
角 的 零 频 点 通过 对 角 变 换 移动 到 频谱 中 间 。 这 里 的 对 角 变换 并 不 是 严格 意义 的 关于 图 片 中 心 的 坐标 及 
灰 度 值 对 调 ， 而 是 将 图 片 的 分 块 区 域 对 调 ， 将 傅 里 叶 变 换 的 图 像 以 图 像 中 心 为 原点 等 分 成 4 AOD, 
外、 加 、 轩 ,通过 将 和 @@、 回 各 进行 区 域 对 调 ， 使 频谱 能 量 移动 到 图 像 中 心 ， 如 图 5(b)。 

将 光斑 质心 连 线 ， 由 连 线 斜率 计算 条 纹 的 角度 。 通 过 连通 区 域 质心 连 线 的 方法 求 出 3 个 光斑 所 在 
直线 的 斜率 ， 求 解 质心 的 方法 是 : 首先 对 图 5(b) 进行 二 值 化 ， 然 后 求 出 每 个 光斑 的 连通 区 域 ,， 计算 
每 个 连通 区 域 的 质心 ， 二 值 化 处 理 时 需要 选择 合适 的 阔 值 。 图 5(c) 是 传 里 叶 对 角 变 换 图 像 的 二 值 化 
图 ， 为 求 连通 区 域 做 准备 。 图 5(d) 是 二 值 化 图 中 的 3 个 连通 区 域 ， 并 用 十 字 丝 标记 每 个 连通 区 域 的 
质心 位 置 ， 这 里 是 以 图 像 左上 和 角 的 第 1 个 点 作为 坐标 原点 。 


根据 斜率 公式 可 得 到 两 质心 之 间 的 斜率 厂 ， 忆 ,三 ， 求 出 质心 连 线 的 斜率 的 平均 值 E= ktkt), 


输出 干涉 条 纹 角度 


图 4 傅 里 叶 变换 计算 条 纹 角 度 算 法 流程 图 
Fig.4 Flowchart of the algorithm for calculating 


fringe angle based on Fourier transform 


所 以 条 纹 的 斜率 = 二， 根据 得 到 的 斜率 计算 出 条 纹 角 度 9。 
k 
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换 ; (b) 对 角 变 换 ; (c) 二 值 化 图 ; (d) 边缘 轮 廊 及 质心 图 


Fig.5 (a) Fourier forward transform; (b) Diagonal transform; (c) Binary image; (d) Contour and centroid map 


图 5 (a) #2; eR 


TERE PAE A A HE, BoA REI, TL BE ACER, seh eR he TI k BE 
的 问题 。 以 8 位 无 符号 数据 为 例 , 在 阔 值 范围 1 ~255 PERAH, PRAIA TENA, i+ 
结果 如 图 6。 对 合理 阔 值 范围 内 的 所 有 角度 求 平 均值 average_Angle = 43. 0078”， 与 理论 值 43 之 间 的 
误差 o=0.018%, 
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Fig.6 Maps for reasonable threshold range of simulation data and experimental data 
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3 实验 干涉 条 纹 算法 验证 


独立 三 孔径 干涉 实验 的 参数 与 仿真 数据 的 参数 基本 一 致 ， 但 因 测量 误 差 和 有 瞳 映射 误差 的 原因 ， 子 孔 
径 透 射 到 合成 镜 上 的 光斑 尺寸 ， 以 及 基线 的 尺寸 与 仿真 数据 存在 一 定 差 异 ， 实 验 采 集 的 条 纹 如 图 7(a) 。 

按照 2. 2 节 中 的 算法 得 到 如 图 6 的 合理 冰 值 范围 ， 对 比 上 下 两 图 可 知 实 验 数据 和 仿真 数据 的 合理 
Bd (HIE A IME SO, 64] ， 在 此 范围 计算 得 到 的 角度 的 平均 值 为 44.790 1 。 从 实验 和 仿真 的 
结果 对 比 来 看 ， 两 者 的 合理 闪 值 范围 并 不 完全 一 致 ， 这 是 由 于 实验 采用 的 光学 需 件 以 及 整个 光学 系统 
会 引入 像 差 ， 以 及 采集 条 纹 的 相机 存在 噪声 ， 而 仿 直 所 用 的 条 纹 图 是 理想 的 无 像 差 、 无 噪声 图 像 。 

通过 步 进 扫描 法 得 到 条 纹 数据 ， 此 时 的 条 纹 有 一 定 角度 ， 为 了 获得 子 孔 径 间 的 相对 光 程 差 的 最 小 
位 置 ， 需 要 沿 条 纹 的 法 线 方向 行 采样 。 在 对 行 采样 数据 处 理 时 ， 条 纹 的 对 比 度 作 为 评价 条 纹 好 坏 的 一 
个 重要 指标 ， 甚 定义 是 


ee ee 
Da + Lmin | 
其 中 ，.,,,、1, 分 别 表示 干涉 条 纹 光 强 的 最 大 值 和 最 小 值 。 为 了 方便 对 比 结果 ， 这 里 给 出 两 种 条 纹 ， 
第 1 种 是 不 经 过 相机 旋转 采集 得 到 的 条 纹 ， 如 图 7(a) ， 根 据 计算 得 到 的 条 纹 角度 44.790 1 ， 沿 条 纹 
的 法 线 方向 (-45. 209 9 ) 行 采样 ， 得 到 对 比 度 的 变换 曲线 如 图 7(b) ， 步 进 扫描 的 对 比 度 的 变化 趋势 
并 不 明显 ， 也 无 法 从 拟 合 曲线 中 找到 最 小 光 程 差 的 位 置 。 


y= (4) 
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图 7 未 旋转 相机 。(a) 实验 采集 条 纹 图 ; (b) 沿 条 纹 法 向 行 采 样 得 到 的 对 比 度 变 化 曲线 


Fig.7 Camera not rotate. (a) Fringe from experiment; (b) The contrast curve along the normal direction of the fringe 


另 一 种 是 相机 旋转 一 定 角度 后 采集 的 条 纹 ， 根 据 计 算得 到 的 条 纹 角 度 44. 790 1”， 先 将 相机 逆 时 
针 旋 转 45. 209 9”， 使 条 纹 竖 直 ( 即 90” ) ， 如 图 8(a) ， 然 后 沿 条 纹 的 法 线 方向 (0 ) 行 采样 ， 得 到 对 比 
度 的 变换 曲线 如 图 8(b)。 

由 图 8(b) 可 知 ， 条纹 从 模糊 到 清晰 的 变化 趋势 比较 明显 ， 且 可 以 从 拟 合 曲线 中 找到 条 纹 最 清晰 
的 位 置 ， 也 就 是 光 程 差 最 小 时 对 应 的 促 动 嚣 位置。 但 从 图 8(a) 可 以 看 出 ， 采 集 的 条 纹 质量 不 高 ， 在 
步 进 扫描 的 过 程 中 ， 两 光斑 在 焦 面 位 置 会 略微 分 开 ， 导 致 看 上 去 有 两 组 条 纹 ， 这 是 由 于 在 步 进 扫描 的 
过 程 中 引入 倾斜 误差 。 

通过 实验 可 知 ， 旋 转 相机 后 更 容易 根据 条 纹 对 比 度 的 最 大 值 找 到 最 小 光 程 差 位 置 ， 这 是 由 采集 条 
纹 图 的 相机 本 身 属性 决定 的 。 相 机 采用 2 048 x 2 048 像素 点 组 成 的 点 阵 ， 条 纹 图 是 用 有 限 个 像素 点 组 
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成 的 点 阵 图 。 理 想 情况 下 ， 如 果 条 纹 水 平 或 竖 直 ， 那么 条 纹 像 素 点 最 大 限度 地 分 布 在 一 条 直线 上 。 当 
条 纹 呈 一 定 角 度 时 ， 描 述 条 纹 的 像素 点 实际 上 是 直线 附近 像素 点 的 近似 取 整 。 因 此 在 计算 未 旋转 相机 
之 前 的 条 纹 对 比 度 时 ， 选 择 的 光 强 最 大 值 与 最 小 值 并 不 精确 ， 从 而 导致 对 比 度 的 计算 不 准确 ， 这 是 旋 
转 相 机 使 条 纹 竖 直 后 可 以 改善 确定 条 纹 对 比 度 位 置 的 物理 原因 。 


4 结 


Contrast 


Actuator Position/um 


(a) (b) 


图 8 相机 旋转 后 。(a) 实验 采集 条 纹 图 ; (b) 沿 条 纹 法 向 行 采 样 得 到 的 对 比 度 变 化 曲线 


Fig. 8 Camera rotate. (a) Fringe from experiment; (b) Contrast curve along the normal direction of the fringe 


论 


通过 仿真 数据 可 知 ， 靖 值 不 同 得 到 的 条 纹 角度 也 有 微小 差异 ， 通 过 浆 值 遍历 计算 平均 值 的 方法 ， 
得 到 的 角度 与 理论 角度 更 加 接近 。 将 此 算法 用 到 三 孔径 干涉 实验 ， 根 据 得 到 的 条 纹 角度 ， 指 导 相 机 的 
旋转 角度 ， 使 条 纹 的 角度 沿 竖 直方 向 ， 便 于 条 纹 行 采样 ， 更 容易 、 更 精确 得 到 对 比 度 最 大 的 位 置 ， 即 


最 小 光 程 差 位 置 。 
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A Method for Calculating the Angle of Interference Fringes Based on 
Fourier Transform and Threshold Traversal of Binary Image 
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Abstract; Eliminating the relative Optical Path Difference (OPD) between two sub-apertures in synthetic 
aperture telescope is a precondition for high resolution image. The fringe detection method is effective to 
decrease the relative OPD. The fringe angle depends on the position distribution of sub-aperture. The precise 
angle of fringe is an essential prerequisite for taking line sample along the normal direction. And then, on the 
basis of the maximum contrast value, the position corresponding to the minimum OPD between two sub- 
aperture can be obtained. Here, a method for calculating the angle of interference fringes based on Fourier 
transform and threshold traversal of image binaryzation is proposed. Firstly, the basic principle of the algorithm 
is introduced. Secondly, with the data for fringe angle of 43 degrees as the simulation data to test this 
algorithm, the result of 43.007 8 degree is obtained. The value of error is 0.018% compared with theoretical 
simulation value. In the end, contrast curves from two cases of rotating camera and non-rotating camera are 
compared. It is more benefit to get the accurate minimum OPD by the method of rotating camera to make the 
fringe along the longitudinal axis of camera according to the angle of fringe. 


Key words; Fourier Transform; Threshold; Interference fringe; Fringe angle; Optical path difference 


